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We solve the quasiclassical problem of tunneling through an external potential barrier in a 
dense plasma, where the tunneling particles undergo simultaneous collisions with other particles in 
thermodynamic equilibrium. Under such conditions the spectral density of states available to the 
particle has a Lorentz shape, rather than the delta function, which leads to a quantum tail in the 
particle momentum distribution function. We show that, this tail indeed significantly alters the 
average nuclear reaction rates, which supports earlier suggestions. This rate can be many orders of 
magnitude higher than would be normally calculated by averaging over the Maxwell distribution of 



energies. 

The particle distribution over momenta, /(p), 
acquires non-Maxwellian tail in thermodynamic 
equilibrium [|IJ due to quantum effects, while the 
particle distribution over energies, /(&>), remains 
Maxwellian. An interesting suggestion was made 
in Ref. 0, that the rates of various processes, 
including nuclear reactions, may be significantly 
increased by the quantum tail of the momentum 
distribution. A final formula for the ionization 
rate in a closed form was obtained in @, by us- 
ing the Born approximation. The formula for 
the nuclear reaction rates was found in || un- 
der assumption that it is correct to substitute 
the particle momentum in the form e p = p 2 /2m 
into the known quasiclassical cross-section a{uS) 
and then average it over the particle distribu- 
tion over momenta /(p) rather than the distri- 
bution over energies f(oo). This procedure be- 
comes even more unclear, if we take into account 
the fact, that the momentum, p, and energy, lu, 
of a colliding particle are independent variables 
not connected by the usual dispersion relation 
5(u — p 2 /2m) |],[|,@. Note that the averaging 
of o~{uj) over the distribution over energies f{w), 
which is Maxwellian would give the usual reac- 
tion rates || rather than the rates found in [[|, 
which are accelerated by many orders of magni- 
tude in certain regimes. It is therefore important 
to find the nuclear reaction rates from first prin- 
ciples. 

In this paper we rigorously solve the quasiclas- 
sical problem of tunneling through the potential 
barrier, when the tunneling particles undergo si- 
multaneous collisions with other particles of the 
plasma maintained in thermodynamic equilib- 
rium, (see Fig. 1). The plasma has density and 



temperature n ,T and is fully ionized. We use 
the Green-function technique and do not 

rely on any assumptions || about any averaging 
procedure. We show that, if the requirements 
of the quasiclassical approximation are fulfilled 
with respect to both the barrier and collisions, 
our final result for the nuclear reaction rate sup- 
ports the idea of averaging over momentum pos- 
tulated in ||. 
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FIG. 1. 



According to the diagram technique of 
Keldsyh and Korenman f7||| the kinetic prop- 
erties of the particles are described by the Green 
function G~ + (X 1 , X 2 ); we adopt the notation 
X = (i, x). Kinetic equations for this func- 
tion J7|,[8| include the interaction between the 
particles as well as their interaction with the ex- 
ternal field. Therefore, the rate of change of this 
function, evaluated at any point after the bar- 
rier, will give the tunneling rate K(x, t) for the 
particles colliding simultaneously with other par- 
ticles: 



K(x,t) = \(d tl -d t2 )G- + (t 1 * 1 -t 2 K 2 ) 



X\ — *X 2 =x 



d tl G + (t 1 x 1 ;t 2 x 2 ) 



Xi^X2=X 



(1) 



In the second line of ([I]) we anticipated that, G~ + 
depends only on t\ — t 2 . 

Now we will proceed with solution of the ki- 
netic equations for G~ + (Xi,X 2 ). They can be 
written as Dyson equations for G 
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where we use subscripts 1..4 to denote depen- 
dence on X1..X4. Dependence on G® 2 can 
be eliminated || by acting on both sides of the 
Eq. (El) with the operator 



G7 1 



2m 



where = + d 2 ± . Since we are interested 
in the steady state solution, one can see that, 
the properties of this operator can be fully ac- 
counted for if we introduce the function g(z, k), 
such that 



2m dz 2 



-U(z) 



g(z,k) = -e k g(z,k), (3) 



where = k 2 /2m. Since the operator L is 
Hermitian, we will use the property of its eigen- 
functions g(z, k): 



g(z,k 1 )g*(z,k 2 )dz = 2ii5(k 1 -k 2 ) (4) 



Now let us solve for Gf 2 + , which satisfies 



v(0) 



0. 



(5) 



It is easier to find by using its definition 



G g)-+ _ , /.r.t. 



(6) 



where ip are the Heisenberg operators of cre- 
ation and annihilation of the particles at a point 
X, and < .. > means quantum and statistical av- 
eraging. Since ip evolves according to G~ l ip = 0, 
one can see, that it is equal to 

${X) = £ a^giz, k )e-^ t+l ^ (7) 

q x fc 



Here a q±/ t is the annihilation operator of the 
particle with momentum (q±,k), and e q± fc = 
q\/2m + k 2 /2m. 

We substitute ip(X) from (0) into Eq. (|6]) and 
obtain 



G(± + ( Zl ,z 2 ) = 2m I ^ 



00 dk 
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6(u - e q±fc + ^)n(q 1 _k)g(z 1 , k)g*(z 2 , k), (8) 

where we denote (a qx fcfl q± fe) = n(q±k); jj is the 
chemical potential. 

It is easy to see that, in thermodynamic equi- 
librium we are considering, the function n(qj_k) 
is equal to Fermi distribution 7ip(e qx fc). To con- 
vince ourselves we take G^~ + (zi, z 2 ) at such 
values of -21,-22, when U{z) is negligible and 
g(z, k) =>■ exp(ikz). One more Fourier transform 
Jexp[-iq(z 1 - z 2 ))d(z 1 - z 2 ) will give G^ + , 
which has to coincide with the known expres- 
sion § G<®-+ = 2m8{u - e q + //)n F (e q ). ' 

Before proceeding further with solving for 
G- + (X U X 2 ) we first find K<®, which is the tun- 
neling rate obtained by using G^ 2 ~ + in Eq. (JI|). 
Since we are interested in behind the bar- 
rier, where U(z) = 0, we can exchange <\ for 
iA\/2m as is clear from the form of the operator 
Gf 1 and Eq. (^). Then the formula for the rate 
takes the form: 



du 
2^ 



d 2 q_i_ 
(2n) 2 



F) 2 

q£2 wqj 



(*,0 



JC=o 
(9) 



where z — {z\ + z 2 )/2, £ — z x — z 2 . 

In deriving (|9|) we expressed G ( f 2 ~ + through 
its Fourier transform G$J.~ + (z, () with respect to 
"fast" space r_|_ = — r% and time r = t\ — t 2 
variables and took the limit r, r — > as pre- 
scribed by Eq. (|1|). 

Finally we need to know the function g(z, k) 
in the region behind the barrier. We will use the 
standard quasiclassical expression M 



g(z, k) 



Gh 



q{z,ky/ 2 



exp 



i / q(z', k)dz' 



(10) 



C k = A k exp 



z b (k) 



q(z', k)dz 



q(z, k) = y/2m(e k -U(z)) 



(12) 



The square of the exponential factor in Eq. ( |TTD 
is the tunneling coefficient, W(k) ||10|| . We multi- 
ply it by a factor A k = (S (e k ) / e k ) 1 '\ which takes 
into account nuclear physics effects, not consid- 
ered here. S(e k ) is the astrophysical factor ||TT |. 

To find from Eq. flj^) we substitute 

Eqs. ©-(0) in Eq. ©, use z x = z + (/2, 
z 2 = z — C/2 and take d 2 derivative in the quasi- 
classical sense keeping only zero-th order terms. 
The result of the differentiation is 



8P_ 



[g(z 1 ,k)g*(z 2 ,k)} 



-q(z,k)C 2 (13) 



J(=o 



Note, that q(z, k) = k is true behind the barrier. 

The problem of tunneling collision between 
particles mi and m 2 can be reduced to the one- 
dimensional motion of a particle with reduced 
mass m r = m\m<il [m\ + m^) in the external po- 
tentianl U (r) || . We make this reduction in the 
original Hamiltonian and then second quantize 
it. Then we use spherical coordinates, assume 
spherical symmetry of the tunneling collision and 
perform calculations similar to those explained 
above. The final result for K^ 1 is 



tAo\ f 00 4nk 2 dk .. . k S(e k ) rir ,,. 
K { ) = K> / -^rr^n M {k) b^W(k), 



o (2vr) 3 



m r e k 



(14) 



where e k = k 2 /2m r . W(k) is obtained from 
W(k) by z —>■ r, n,M{k) is the Maxwell distri- 
bution which depends on m r . 

It agrees with the answer obtained in || 
by averaging the quasiclassical tunneling fac- 
tor W(k) over the Maxwell distribution, 
{riMiy )a{m r v 2 /2)v) . Hence, we see, that the dis- 
persion relation 5{uo — e q ) was implicitly used in 
||, which corresponds to an approximation of 
instantaneous, two-body collisions. Now we will 
not make this assumption and proceed with find- 
ing the tunneling rate, K, of the particles, which 
always collide and hence are never "in between 
collisions" . 



By using ||£|, we find 



(15) 



where we denote T, R — T, A 
«7o;qx(?> z )- Note, that 7ujq(-2) is a non-linear 
function of G af3 [[§§. We now analyze Eq. (|TJ) 
qualitatively, which will help us to find that part 
of the solution, which makes the largest contri- 
bution to K. 

In the region away from the barrier we can ne- 
glect the LHS of Eq. ([15]), and obtain: 



= -n F {u) (G^(z)-GX( 



(16) 



Therefore, Eq. (|T5|) will allow us to propagate 
this solution into the region behind the barrier. 

We can find the first integral and formally inte- 
grate Eq. ( |1~5D along trajectories, which will lead 
to a sum of solutions of homogeneous and "in- 
homogeneous" equations (due to the last term 
in (jlBD ). As can be seen from Eq. (|i"B"l), "inho- 
mogeneous" solution will involve the integral of 
the product of 7^(2) and G R (z) or G R * = G A . 
Qualitatively, it means, that in the region of 
z behind the barrier it will be proportional to 
a product of at least two exponential factors 
W{k). This is so, because both 7o;q(z) and 
G^ q (z) will depend on the integral of the prod- 
uct g(zi, k)g*(z2, k) and such product is oc W{k) 
as can be seen from Eqs. (p~Q|) - (P~2|) - 

It is obvious that, the solution of the "homoge- 
neous" equation, Gj + , involves only one factor 



of W(Jfe), since Gj + oc gg*, see Eq. (|T3|) . There- 
fore, we will find only the solution to the "homo- 
geneous" equation, which we will use to obtain 
K. We supplement this homogeneous equation 
with the boundary condition, G h from Eq. ( |T5] ) 
imposed at z = z* away from the barrier, see 
Fig. 1. 

Since the largest contribution to K is made 
by Gj + one can see, that we can use Eq. (|9|) to 
find K if we substitute Gj + instead of G^°^ + 
in (||D. Therefore, we find Gj + with dependence 



on Zi, z 2 from the very begining by applying the 
transformation 

dki dk 2 



g(z 1 ,k 1 )g*(z 2 ,k 2 



2tt 2tt 



to Eq. (Q). We use the boundary condition, 
Eq. ( |i~6|) and obtain the following answer: 

/dk 
g k (zi)6^(u - e n±k )g* k (z 2 ) — 

(18) 

where g k (z) = g(z,k) and we used the well- 
known equilibrium solution [Q 

n R (k) - Cr A (k) - 7^ q± ( fc ) 



(19) 



(k) = 5^{uj - e qxfc ) 

M), (PI) we see, 



Comparing Eq. (^|) with 
that the effect of collisions, which makes the 
largest contribution to K can be described as 
8{uj— e q± fc) =>- <5 7 (c<j— e q± fc), as can be expected on 
the intuitive grounds. Note also, that we can ob- 
tain the answer fll8|) , (|l9|) by following the same 
steps leading to Eq. flD if we make a substitution 
e q ± fe — *■ eq ± fc+^7 q± fc- This means that, we have to 
use damped oscillators instead of undamped ones 
for second quantization of colliding particles. 

Now we substitute Eqs. ([18]), ([ID into © and 
integrate / doj/2^ by using the result of [l|]: 



^n M (w)i5 7 (w 

Z7T 



e q± fc) = w M (e q±fc ) + 5n 7 (q±A;), 

(20) 



where <5n 7 (q_i_/c) is the "quantum tail". 

At this point we are ready to perform the same 
steps as those leading to Eq. fll4[ ) in order to de- 
scribe the tunneling of particles mi and m 2 . The 
answer for K is: 



K = n 



A-np 2 dp S(e p ) 
(2tt)3 e p 



V 



W{p) [n M (p)+6n, 



K M + K T (21) 



We now evaluate S R through collision fre- 
quency, ut, as in [0,0, carry out the integral 
in (|2"T[ ) and form the ratio 7*12 = K^/Km'- 



Tl2 



3 19 / 2 ^ hu T (m coll ) f m coll y/ 2 e™ 



T 8 



(17) Here m coH 



m r mj/(m r + m 



1/3 



m im2 /( mi +m 2 ), r 12 = 3(tt/2) 2 / 3 ( 100 ^f Al2 
Ai 2 = A 1 A 2 /(A 1 + A 2 ), T k is the temperature in 
keV; Ax, A 2 are the atomic numbers of tunnel- 
ing particles, irtj is the mass of the background 
particles colliding with tunneling particles. Sum- 
mation is over all species of colliding background 
particles, including tunneling species. 

This result, (PH), (|2"2"|), supports the idea of av- 
eraging of (TV over momentum distribution pos- 
tulated in ||, but differs from the answer of 
by a factor (m co u/m r ) 3 under the summation 
sign. 

Consider, for example, a DT plasma of density 
p = 10g/cm 2 and temperature T = 0.1 keV. 
Then Eq. (^) will lead to the rate acceler- 
ated by five orders of magnitude as compared 
to the conventional answer. We find the ra- 
tio r DT = Kl T /K^ T = 8 • 10 4 , with K% T = 
2.7 • 10- 30 cm 3 /sec, K^ DT = 2.1 • 10- 25 cm 3 /sec. 
In case of heavier elements and not very small 
hPr/T, ratio becomes much higher ||. 
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